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I. Introduction 



The F + H 2 reaction (and its isotopomers F + D 2 and F + HD) has served as one of most 
important benchmark elementary chemical reactions for at least the last four decades [1-4]. 
During this period, the F + H 2 (D 2 , HD) interactions have been the subject of extensive studies, 
both experimental, mainly in crossed molecular beams, and theoretical. The first molecular 
beam data on the dynamics of these reactions was published by Lee and co-workers in 1970 [5]. 
In the mid-eighties the same group reported the vibrationally resolved angular distributions 
of the F + H 2 (D 2 , HD) reaction products at various collision energies E co \ [6-8]. 

In their milestone paper [7], Neumark et al. measured the vibrationally resolved center- 
of-mass (CM) differential cross sections (DCSs) of the F + H 2 (t> = 0; j = 0, 1,2) reaction at 
collision energies of E co \ = 1.84, 2.74, and 3.42 kcal/mol. The most intriguing observation of 
Lee and co-workers [6, 7] was an unexpected noticeable forward peak in the angular distribution 
of the HF(w' = 3) product [9], which increases in intensity with increasing collision energy. 
The angular distributions of the HF(t/ = 2) and HF(t>' = 1) molecules do not exhibit such 
a peak. A similar but less pronounced peak can be seen in the CM DCSs of the DF(t>' = 4) 
product from the F + D 2 reaction [8]. The presence of these peaks has been confirmed in 
subsequent crossed beam experiments for the F + H 2 reaction [10] as well as for the F + D 2 
reaction [11, 12]. 

The vibrationally selective forward peaks in the angular distributions of the HF(t>' = 3) 
and DF(t> ' = 4) products from the F + H 2 and F + D 2 reactions, respectively, were attributed 
by Neumark et al. [6-8] to quantum mechanical (QM) Feshbach resonances, i.e., metastable 
states formed on the vibrational adiabatic potentials in the potential energy surface (PES) 
transition region [4, 13-15]. It is interesting to note that some resonances in the F + H 2 and 
F + D 2 reactions were first theoretically predicted three and a half decades ago [16] in collinear 
collisions. However, subsequent to their observation these peaks for both the reactions were 
reproduced in quasiclassical trajectory calculations on various PESs of the 1 2 A' ground state 
of the FH 2 system [12,17-21]. In view of the classical nature of these calculations this was 
regarded as evidence against an explanation of the peaks by a quantum resonance [3,22]. 
Moreover, Castillo et al. [23] concluded, from an analysis of some QM characteristics of the 
F+H 2 scattering, that the forward peak in the angular distribution of the HF(t> ' = 3) molecules 
results from tunneling through the combined centrifugal and potential energy barrier at large 
values of the total angular momentum J rather than from a resonance. On the other hand, 
the data of several subsequent studies, both experimental [24] and theoretical [25-27], favored 
the original resonance explanation proposed by Lee and co-workers [6,7]. This controversy 
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has also been discussed in several reviews [15,28]. Recently, a new crossed molecular beam 
experiment on the F( 2 P 3 / 2 ) + H 2 (i> — j = 0) reactive scattering was carried out by Wang 
et al. [10] in the collision energy E co \ range from 0.4 to 1.2 kcal/mol. Based on a detailed 
QM simulation they concluded [10] that the forward peak in the angular distribution of the 
HF(V = 3) product is generated by a slowing-down during passage over the centrifugal barrier 
in the exit valley, with a small contribution from a shape resonance [4, 14] at i? co i slightly 
above 0.5 kcal/mol, and, moreover, that Feshbach resonances do not contribute to the forward 
peak. 

In crossed beam experiments of 2006, intense forward peaks were discovered in the angular 
distributions of the HF(t> ' = 2) product from the F( 2 P 3 / 2 )+H 2 (t> = j = 0) reaction at a collision 
energy of E co \ = 0.52 kcal/mol [29] and the F( 2 P 3 / 2 ) + H 2 (t> = 0; j = 1) reaction at a collision 
energy of E co \ = 0.19 kcal/mol [30], much lower than sampled by Neumark et al. in the eighties 
[6, 7]. In contrast to the forward peak of the HF(V = 3) molecules discussed above, the peaks 
pertaining to HF(V = 2) are unanimously regarded as consequences of Feshbach resonances 
[4,29-33]. Moreover, in both the reactions F + R 2 (v = j = 0) and F + H 2 (v = 0; j = 1), the 
forward peak in the angular distribution of the HF(V = 2) product is caused by two dynamical 
resonances which interfere with each other. 

Numerous convincing manifestations of QM Feshbach resonances have also been found at 
integral and differential cross sections of the F + HD — > D + HF reaction at various collision 
energies. The resonances in this reaction have been extensively studied both experimentally 
and theoretically since 2000 [24, 34], as discussed in the reviews [4, 15, 28, 32, 35] and references 
therein. In the last four years several additional important publications have appeared [36-39]. 
The F + HD interaction is interesting from the viewpoint of its stereodynamics [40,41]. 

In this article, we reexamine the controversial role of resonances in the forward scattering 
of the HF(V = 3) product from the F + H 2 reaction employing an approach first described 
in our preceding publications in Russian [42]. Our QM simulations of the F + H 2 interaction 
at the collision energies of the experiment by Lee and co-workers [7] (E co \ = 1.84, 2.74, and 
3.42 kcal/mol) explore the trends in the behavior of the vibrationally and rotationally resolved 
DCSs of the HF(V, j') scattering as the vibrational quantum number v' increases from to 3. 
Although the cross sections for the formation of the HF product from the F + H 2 reaction in the 
ground vibrational state are very small, HF(t> ' = 0) molecules have been detected in a crossed 
molecular beam experiment [43]. Our analysis shows that the vibrationally specific forward 
peak of the HF(V = 3) product can be explained by the superposition of two independent 
effects which reinforce each other. One of these effects is a purely energy restriction, while the 
other one touches upon all the vibrational states v' of HF molecules and is therefore hardly 
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of mainly resonance nature. Thus, this study provides additional evidence in support of the 
conclusions of Wang et al. [10] against a resonance origin of the peak in question. 

The paper is organized as follows. After a brief discussion of the calculation details in 
section II, we analyze the QM vibrationally and rotationally resolved DCSs of the F + H 2 (t> = 
0; j = 0, 1, 2) reaction in section III. The remarks of section IV conclude the article. 

II. DCS Calculations 

The ab initio Stark-Werner (SW) PES [44] of the 1 2 A' ground state of the FH 2 system 
has had an enormous impact on the theoretical development since 1993 [3,45]. Considerable 
effort has since been devoted to further improving this surface in the entrance valley which 
has led to the HSW PES [20,46], SW-LR PES [47], SW-LR-SO PES [47], and PES III [48], or 
in the exit valley (the SWMHS PES [49]), or in both the valleys simultaneously (the PES IV 
[38]). In addition, two totally new ab initio surfaces XXZ [29,50] and FXZ [51] were recently 
produced. The calculations of the present paper are carried out on the standard SW PES 
[44] to facilitate comparison with the data of the previous publications [3,20,23,26,47,52]. 
The SW surface is still widely used in studies of the F + H 2 (D 2 , HD) reactions [40,53-55]. 
Moreover, our observations are of a qualitative character, and we have not attempted a direct 
comparison with the experiment. Recall, however, that the experimental forward peak in the 
angular distributions of the HF(t> ' = 3) product from the F + H 2 reaction is less pronounced 
than that obtained in QM simulations on the SW PES [3,20,23,47,52]. The same situation 
holds for the forward peak in the angular distributions of the DF(i/ = 4) product from the 
F + D 2 reaction [21]. 

As in most of the previous theories of the F + H 2 (D 2 , HD) reactive scattering [3], contri- 
butions from the excited state F*( 2 Pi/ 2 ) of the fluorine atom reactant were neglected. This 
approximation is fully justified for collision energies E co \ > 1.84 kcal/mol in the context of the 
present study [36,52,56,57]. Also according to Tzeng and Alexander [52], these contributions 
affect mainly backward scattering. Note, however, that several ab initio multi-state surfaces 
for all of the three lowest electronic states 1 2 A', 2 A", and 2 2 A' of the FH 2 system as well as for 
the 1 2 A' <H- 2 2 A' nonadiabatic coupling have been constructed: the ASW [56], MASW [58], 
LWA-5 [36], and LWA-78 [36] PESs. 

For each collision energy E co \ = 1.84, 2.74, and 3.42 kcal/mol [7], the vibrationally and 
rotationally resolved DCSs da v 'j//dQ of the F + H 2 (t> = 0;j) — > H + HF(i/, j') reactions were 
calculated for j = 0, 1, and 2 using the ABC program [59]. This code solves the Schrodinger 
equation for the motion of the three nuclei on a given PES by a coupled-channel method 
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in the Delves hyperspherical coordinates [60]. The calculations were performed for the total 
angular momentum J ranging from to a maximum value J max . For i? co i = 1.84 kcal/mol, 
the convergence parameters of the ABC program were set at J max = 25, E max = 1.7 eV, 
jmax = 17, k max = 4, p max = 12 bohr, and M tr = 150, and for E co \ = 2.74 and 3.42 kcal/mol 
at J max = 30, £ max = 2.5 eV, j max = 21, k max = 5, p max = 12 bohr, and M tr = 200. Here 
E max is the maximum internal energy of the H 2 reagent and HF product admissible in the 
basis functions, j max is the maximum rotational quantum number of the H 2 reagent and HF 
product, A; max is the maximum (in absolute value) helicity quantum number of the H 2 reagent 
and HF product, p max is the maximum hyperradius p of the system used while solving the 
hyperradial coupled-channel equations, and M tr is the number of propagation sectors involved 
in solving those equations (for details, see the paper by Skouteris et al. [59]). Some test 
calculations with other values of the convergence parameters indicate that the chosen values 
of the parameters yield the DCSs for almost all the BF(v',j') states with a relative accuracy 
of better than ~ 1%. 

III. Results and Discussion 

a. Forward Scattering Coefficients 

As an example, Figure 1 presents the rotationally unresolved DCSs dcv/df2 and some 
rotationally resolved DCSs davj' /dQ of the F + H 2 (t> = 0; j = 2) reaction at a collision energy 
of E co \ = 3.42 kcal/mol. As a whole, the DCSs depend on E co \ only rather weakly. There 
is, however, a distinct shift towards the smaller 9 values as the collision energy increases. 
In the second panel from the bottom of Figure 1, one observes a moderate forward peak in 
the angular distribution of the HF(t>' = 2) molecules. On the other hand, the peak in the 
angular distribution of the HF(i/ = 3) molecules in the region of small 9 < 20° is very strongly 
pronounced in the lowest panel of Figure 1. For other values of collision energy E C(A and initial 
rotational quantum numbers j of the target H 2 molecule, this peak is also largely confined to 
the CM scattering angle range 0° < 9 < 20° in agreement with the results of the previous 
QM simulations on the SW PES [3,20,23,47,52]. Consequently, to characterize the relative 
contributions from various j' values to the forward peak, each (v',j') state of the HF product 
from the F + H 2 (t> = 0; j) reaction will be assigned the quantity 




(1) 
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which we call the forward scattering coefficient of the HF(V, j') molecules. This coefficient can 
be either positive or negative. The larger the positive I v >j> coefficient, the more pronounced 
the forward peak in the HF(V, j') angular distribution. However, since in many cases the 
rotationally resolved DCSs of the F + H 2 (t> = 0; j) reaction oscillate as 9 — > (see Figure 1), it 
is more appropriate to speak here of an "average" forward angular distribution. The forward 
scattering coefficients of eq 1 involve the integrals / (do~ v ' y j dO) (9) d9 instead of the expres- 
sions f(da v 'j'/dfl) (9) sin 9 d9 since otherwise the DCS behavior at small angles 9 would be 
suppressed. The integrals in eq 1 just compare the mean values of the DCSs over the corre- 
sponding angular ranges and do not measure the scattering into these angular ranges. The 
values of the l v iy coefficients for various E co \, j, v', and j' are presented in Figure 2. 

For v' = 0, all the I j> coefficients are tiny for all E co \ and j except for I 0j/ with j' > 14 
at E co \ = 3.42 kcal/mol and j = 1. The same is also essentially true for v' = 1, but the hj> 
coefficients for j' < 5, although very small, are noticeably larger on the whole than the I\y 
coefficients for j' > 6. For v ' = 2, many of the I 2 j> coefficients for f < 5 are already rather large, 
whereas most of the I 2 j> coefficients for j' > 6 remain small. Finally, for v' = 3, conservation 
of energy restricts the rotational levels of the HF(V = 3) molecules. The corresponding I 3 j> 
coefficients for j' < 5 are on the whole much larger than the I 2 j> coefficients. Note that the 
smaller j and the higher E co \ are, the larger are the l3j> coefficients for any fixed j' < 2. The 
same trends persist, in general, for j' > 3. For this reason, lower initial rotational excitations 
of the target H2 molecule and higher collision energies lead to a more strongly pronounced 
forward peak of da^/dQ. 

The forward scattering coefficients I v >j> compiled in Figure 2 provide the following explana- 
tion for the origin of the forward peak of the HF(t> ' = 3) product. As the vibrational quantum 
number of the HF molecules increases from v' — to v' — 3, the I v >j> coefficients grow rapidly 
for small j' (j f < 4 or j' < 5), whereas for larger j' they do not tend to increase and remain 
small for all v'. The fact that the forward scattering peak of the dcv/dfi DCSs is present 
for v' — 3 and absent for v' < 2 is caused by the joint influence of two effects which are 
independent but nonetheless strengthen each other. 

First we note that among the rotationally resolved DCSs da v ff/dQ, only the cross sections 
corresponding to small j' values possess such a peak. However, for v' < 2, the HF(v',j') 
products with small j' constitute only a fraction of all the HF(V) molecules in the given 
vibrational state v' [61]. Consequently, the presence of the forward peak of the dcr v >j>/dQ 
DCSs with small j' values cannot substantially affect the total angular distribution of the 
HF(V) molecules. For v' = 3, on the other hand, as pointed out above, all the BF(v',j') 
products have small rotational quantum numbers j' due to energy restrictions and exhibit 
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forward scattering peaks. After the summation over j', these peaks yield a forward peak in 
the rotationally unresolved angular distribution of the HF(V = 3) molecules. 

Secondly, as v' increases, the stronger are the forward peaks of the rotationally resolved 
DCSs da v iji/dfl for small f. For this reason, even the sum 

3 <5 

with all the rotational levels j' > 6 excluded, exhibits for v' = 3 a much more pronounced 
forward peak than for v' = 2, not to mention v' < 1. 

Thus, the forward peak of the HF(t> ' = 3) molecules discovered experimentally by Neumark 
et al. [6, 7] is most likely explained by the different behaviors of the HF DCSs not only with 
respect to the vibrational quantum number v' but also with respect to the rotational quantum 
number j'. The first of the effects indicated above, the absence of the HF(V = 3; j') products 
with large j', is just an energy restriction. The second effect, an increase in the forward 
scattering peaks of the rotationally resolved DCSs da v >j//d£l for small j' as v 1 grows, is a trend 
that affects all the vibrational states v' of the HF product. The resonance origin of this effect 
seems therefore rather questionable. We conclude that the analysis of the forward scattering 
coefficients I v /ji provides evidence against the resonance nature of the forward peak of the 
HF(V = 3) molecules. 

For the rotationally resolved DCSs da v >j> / dQ the first effect is of no consequence compared 
to the rotationally unresolved DCSs da v >/dQ. The increase in forward scattering of the rota- 
tionally resolved HF products in passing from v ' — 2 to v' — 3 is therefore much less sharp 
than that of the rotationally unresolved products as illustrated in Figure 1. The rotationally 
unresolved DCSs da v > /dO (solid black curves in Figure 1) have a forward peak for the da^/dfl 
cross section and with the exception of a small rise in the da 2 /dQ cross section no such peak for 
v' < 2. On the other hand, the rotationally resolved DCSs dcv^/dfi show a distinct gradual 
increase in the forward peak for small j' values in passing from v' — to v' — 1, v' — 2, and, 
finally, to v' — 3. 

A similar situation holds for the F + D 2 reaction. Figures 2-4 in the paper by Martmez- 
Haya et al. [21] present some rotationally resolved DCSs dcr v >ji/dQ of the F + D 2 (i> = 0; j = 
0,1,2) reaction for v' = 2, 3, and 4 at the collision energies of the Gottingen experiments 
[12,62] E co \ = 2.08, 3.23, 4.15, and 5.53 kcal/mol. These DCSs, averaged over j according 
to the experimental rotational distributions of the target D 2 molecules, were obtained from 
quasiclassical trajectory as well as from quantum mechanical close-coupling calculations, both 
on the SW PES. For the QM cross sections at E co \ > 3.23 kcal/mol and for the quasiclassical 
cross sections at E co y > 4.15 kcal/mol one finds the same trend as in Figure 1 of the present 



7 



work, i.e., a relative increase in the forward peaks of the d<7</j//df2 DCSs in passing not only 
from v' — 3 to v' — 4 but also from v' — 2 to v' — 3. Unfortunately, Martinez- Haya et al. [21] 
did not pay due attention to this phenomenon. The increase in forward peaks in the angular 
distributions of the DF(v',j') molecules in passing from v' — 3 to v' — 4 is, however, much 
sharper than in the angular distributions of the HF(v',f) molecules in passing from v' — 2 to 
v' = 3. 

The gradual evolution of the forward scattering coefficients as v' increases is enhanced if the 
definition of these coefficients is changed so that the backward scattering of the HF products, 
which is very strong for v' < 2 (see Figure 1), is removed. This is achieved by replacing the 
integration over d9 from 0° to 180° in the denominator of the right-hand side of eq 1 by an 
integration from 0° to 90°: 

r 20° A a , ., r 40° A a , ., 
rnew _ 1nn Jp° (111 J20° QiZ /n\ 

'v? - iUU ,90° ru ; ., • y z ) 

L ^ M 

These modified forward scattering coefficients are presented in Figure 3. Their behavior differs 
but slightly, as a whole, from the behavior of the coefficients of eq 1, but the absolute values 
are considerably larger for v' < 2. Compared with the I v /ji coefficients, the increase in the 
coefficients of eq 2 for small j' levels in passing from v' — 1 to v' — 2 and from v' — 2 to v' — 3 
is much weaker. 

Note that the forward peaks of the experimental DCSs for the HF molecules with the 
maximum vibrational quantum number v' — 3 are confined to the interval 0° < 9 < 30° [7]. 

b. Contributions from Separate Rotational States 

To gain further insight into the forward peak in the angular distribution of the HF(t> ' = 3) 
molecules, we also define the quantities 

r 20° pin- , ■, 

Dv, = 100^^ , (3) 



/o° dn 

so that 

Dy'jl = 100 

j' 

for any v' . The D v <y ratios characterize explicitly the percentage contributions of the f 
rotational levels to scattering of the HF(f') product into the angular range 9 < 20°. The 
values of the quantities of eq 3 for various E co \, j, v', and j' are presented in Figure 4. 
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As is seen in Figure 4, the forward scattering of the HF(i>') molecules is either rotationally 
hot or "neutral" for v' — or rotationally cold for v' > 1. The larger the v' is, the smaller are 
the rotational quantum numbers f of the HF molecules which make a major contribution to 
the HF(f') forward scattering. For v' > 2, this trend still enhances the role of the HF(V, j') 
molecules with j' < 5. The dominant contribution to the forward peak of the da^/dfl DCSs 
comes from the HF(V = 3; f) products with very small f values, namely, with f < 3. The 
only exception is the F + H 2 (i> = 0; j = 2) reaction at collision energies E co \ > 2.7 '4 kcal/mol, 
for which the j' = 4 must also be accounted for. This result is in agreement with the data 
from previous simulations [18,20]. Note that the forward peak in the angular distribution of 
the DF(V = 4) product from the F + D 2 reaction is also rotationally cold [12, 21]. 

The calculations show that the quantities of eq 3 would not be significantly affected on the 
whole if the factor sin 9 is removed from the integrand in both the numerator and denominator. 

c. DCS Oscillations 

The oscillations of the DCSs of the HF(v',j') forward scattering are clearly seen in Figure 1. 
They are especially pronounced for j' = 0, and for v' > 2 their amplitudes increase sharply 
as 9 — > 0. The angular "periods" of the oscillations lie between 8° and 20° in typical cases. 
These oscillations can be attributed to the oscillations at large values of the total angular 
momentum J of the reduced d{, k entries of the Wigner rotation matrix [14,63]. Here k is 
the helicity quantum number of the H 2 reagent while k' is the helicity quantum number of 
the HF product [14,64]. This oscillatory structure is only present in forward and sometimes 
sideways HF scattering. In backward scattering the major contribution comes from small J 
values for which the d{, k (ir — 9) are non-oscillatory functions. Oscillations of the rotationally 
resolved angular distributions of the BF(v',j') products from the F + H 2 reaction have been 
reported previously in theoretical papers [10,20,33,42] but have as far as we are aware not 
been observed experimentally. 

The oscillatory structure can be described quantitatively as follows. Let Mi, M 2 , . . . denote 
the values of the successive maxima of the function (da v ij> /dil)(6) as 9 increases from 0° to 40° 
and TOi,m 2 , . . . the values of the successive minima. Forward scattering at 9 = 0° is regarded 
as a maximum point if 

and as a minimum point if the opposite inequality 
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holds. In the first case (da v /j>/dQ)(0°) = M\ and in the second case (da v >j> /dQ)(0°) = mi. 
The choice of 0.5° in the above inequalities is arbitrary and any other very small positive angle 
could be used. The quantity 

= 10Q max{M 1; M 2 }- minimum,} 

dcvj/ 
max J (9) 

0°<6K90° d\L 

will be called the relative oscillation amplitude for the DCSs of the HF(V,j') forward scat- 
tering. If the da v ij>/dQ cross section has less than two maxima or less than two minima in 
the interval 0° < 9 < 40°, oscillations are assumed to be absent and we set A v >ji = 0. The 
values of the relative oscillation amplitudes A v >j> for various E co \, j, v', and j' are presented in 
Figure 5. 

Figure 5 shows a distinct increase in A v iy for small j' < 5 as v ' grows from to 3. Moreover, 
the relative amplitudes A v /j> for v ' = 1 and v' = 2 decrease in general as j' increases, and the 
values of these amplitudes for the F + H 2 (i> = 0; j = 2) reaction are smaller, as a whole, than 
those for the F + H 2 (t> = 0; j) reactions with j = and j = 1. The reason is that the HF(V, j') 
products are formed only in collisions with helicity quantum numbers \k\ < j and \k'\ < j'. 
Since the oscillations of d{, k functions with different k and kl damp each other, the oscillations 
of the HF(V, j') angular distributions are in general more pronounced for the smaller j and 
j'. Why some da v >ji/dQ DCSs with j' = 1 exhibit no oscillations is at present not clear. 

If the maximum over the interval 0° < 9 < 90° in the denominator of the right-hand side 
of eq 4 is replaced with the maximum over the interval 0° < 9 < 180°, the general behavior of 
the relative oscillation amplitudes would change only slightly, similarly as in the case of the 
forward scattering coefficients (eqs 1 and 2). However, for v ' < 2 and especially for v' < 1, 
this replacement considerably reduces many of the relative amplitudes. 

d. Partial Wave Analysis 

To examine one more facet of the HF(V, j') angular distributions, we introduce another 
quantity denoted by X^j,, which is the partial DCS of the HF(V, j') product scattering (at 
fixed E co \ and j) obtained by taking into account the total angular momenta from to some 
value J < J max - Since the da v iji /dVt cross sections are calculated using the total angular 
momenta from to J max , we have X^r = da v /ji/dil. The contribution to the (do\,y /dfi) (9) 
DCS from a single partial wave corresponding to a given value J of the total angular momentum 
can then be defined as the difference X^,{9) — X^,~}{9) [10,23,53], where for J = we set 
X~,j, = 0. Because of interference effects, such a contribution can be negative for some angles 
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9. The total contribution from the J partial wave to the HF(v',j') forward scattering can be 
measured by the integral 

For any set of values of E coh j, v', and /, a small fraction of negative Yj jt integrals among all 
the nonzero numbers Y®,j,, Y^i, . . . indicates mainly constructive interference of the partial 
waves in forward scattering of the EF(v',f) product. In the region 9 < 20°, almost every 
partial wave reinforces the sum of waves corresponding to smaller momenta J. Conversely, 
a large fraction means mainly destructive interference. In the region 9 < 20°, many of the 
partial waves attenuate the sums of waves corresponding to the previous values of the total 
angular momentum J. 

Our calculations show that for v' — 0, the fraction of negative Yj jf quantities lies between 
20% and 60% for most of the values of E co \, j, and j'. For v' — 1 and v' = 2, it typically lies 
between 20% and 50% and between 10% and 45%, respectively. For v' = 3, this fraction does 
not exceed 35% in all the cases and does not exceed 25% for j' < 4. Moreover, for many of 
the sets (E co \, j, j'), all the integrals Y^, are non-negative. The total fraction of negative Yjj, 
quantities over all the values of E co \, j, j', and J is equal to 40.4%, 36.7%, 27.1%, and 7.2% for 
v 1 = 0, 1, 2, and 3, respectively. Thus, partial wave interference in HF(V) forward scattering 
is increasingly constructive for the larger v 1 (for details, see the Russian paper by Azriel et 
al. [53]). An almost completely constructive character of the partial wave interference for the 
HF(V = 3) scattering in the forward direction was first found by Castillo et al. [23], but their 
analysis was limited to the rotationally unresolved angular distributions and was not carried 
over to values of v' < 2 [23]. 

IV. Conclusions 

As we have seen, many features of the rotationally resolved DCSs da v 'j> / df2 for the HF(V, j') 
product from the F + H 2 reaction change monotonically as the vibrational quantum number 
v' grows from to 3 (and it is one of the main goals of the paper to attract attention to 
this phenomenon). For instance, the forward scattering peak increases (section Ilia), the 
forward scattering cools down rotationally (section Illb), the QM oscillations of the angular 
distributions become more pronounced (section IIIc), and the partial wave interference in the 
small 9 region becomes increasingly constructive (section Hid). Undoubtedly, these trends also 
hold for the F + D 2 reaction. Since they affect all the vibrational states of the HF product 
they seem to hardly arise from a QM resonance. On the other hand, it is the combination of 
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such effects and the energy limitation on the formation of HF(V = 3) molecules with large 
f values that favors the forward scattering peak of the HF(V = 3) products. Thus, most 
probably, QM resonances do not play a key role in the origin of this peak, in agreement with 
some previous papers [10,23]. 

This conclusion (arrived at without handling resonances themselves) is confirmed by some 
other facts. The forward scattering peak is observed in the vibrationally resolved angular 
distributions of the product molecules with the maximal possible v' value from both the F + H 2 
reaction and the isotopically substituted F + D 2 reaction [7,8]. It exists in a rather wide range 
of the collision energies E co \ for any value j = 0, 1, and 2 of the rotational quantum number 
of the diatomic reactant. Moreover, as was mentioned in section I, this vibrationally selective 
forward peak can be reproduced for both the reactions in quasiclassical trajectory calculations 
[12, 17-21]. Note that trajectory simulations fail to reproduce the resonance patterns of the 
F + HD -> D + HF reaction [15, 24, 32, 34, 35]. 

We believe that our approach based on a careful examination of the state-to-state DCSs for 
all the values of v' and j' could be useful in other situations as well, in particular, in studies 
of indisputable resonance effects [4,24,29,30,34,35]. 

It is interesting to note that in the very recent paper by Xiahou and Connor [54] (which 
also uses the SW surface), the helicity-resolved DCSs of the F + H 2 (t> = j — 0) — >■ H + HF(t>' = 
j' = 3) reaction at a collision energy of 2.74 kcal/mol are found to be an example of broad 
(attractive) rainbow scattering. 

A remaining important task is to clarify the relationship between the features of the PES 
topography and the dynamical characteristics of the F + H 2 reactive scattering discussed in this 
work, for instance, the dominance of the positive forward scattering coefficients I v ij> over the 
negative ones and their growth as v' increases for fixed small j' values [65, 66]. The correlation 
techniques we introduced recently [66, 67] could probably be of help here. However, the relative 
trends in the DCS behavior studied in this work cannot be highly sensitive to the details of the 
surface used. In particular, all these trends are expected to also hold for the more advanced 
ground state FH 2 PESs [20,29,38,46-51] listed in section II. 
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Figure Captions 



Figure 1. Some vibrationally and rotationally resolved DCSs da v >j'/dQ of the F + H 2 (i> = 
0; j = 2) reaction at the maximum collision energy E co \ = 3.42 kcal/mol. 

Figure 2. The forward scattering coefficients l v iy of eq 1. 

Figure 3. The modified forward scattering coefficients I^fJ, of eq 2. 

Figure 4. The contributions D v >j> of eq 3 from the rotational levels j' to forward scattering 
of the HF(V) products. The meaning of the curves is the same as in the v' = 1 panels of 
Figures 2 and 3. 

Figure 5. The relative oscillation amplitudes A v iy of eq 4. The meaning of the curves is the 
same as in the v' — 1 panels of Figures 2 and 3. 
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